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Abstract. About one-quarter of the extrasolar giant 
planets discovered so far have orbital distances smaller 
than 0.1 AU. These "51 Pegb-like" planets can now be di- 
rectly characterized, as shown by the planet transiting in 
front the star HD209458. We review the processes that 
affect their evolution. 

We apply our work to the case of HD209458b, whose 
radius has been recently measured. We argue that its ra- 
dius can be reproduced only when the deep atmosphere is 
assumed to be unrealistically hot. When using more realis- 
tic atmospheric temperatures, an energy source appears to 
be missing in order to explain HD209458b's large size. The 
most likely source of energy available is not in the planet's 
spin or orbit, but in the intense radiation received from 
the parent star. We show that the radius of HD209458b 
can be reproduced if a small fraction (~ 1%) of the stellar 
flux is transformed into kinetic energy in the planetary at- 
mosphere and subsequently converted to thermal energy 
by dynamical processes at pressures of tens of bars. 

Key words: extrasolar planets, 51 Peg, HD209458, evo- 
lution, interior structure, giant planets 



1. Introduction 

The detection of planetary-mass companions in small or- 
bits around solar-type stars has been a major discovery of 
the past decade. To date, 73 extrasolar giant planets (with 
masses Msmi < 13 Mj, Mj being the mass of Jupiter 
and i the inclination of the system) have been detected by 
radial velocimetry. Fifteen of these (21%) have distances 
less than 0.1 AU, and ten (14%) have distances less than 
0.06 AU (see Marcy et al. |2000| and the discoverers' web 
pages) . This is for example the case with the first extraso- 
lar giant planet to have been discovered, 51 Peg b (Mayor 
& Queloz 1995). These close-in planets form a statisti- 
cally distinct population: all planets with semi-major axis 
smaller than 0.06 AU have near-circular orbits while the 
mean eccentricity of the global population is < e >« 0.27. 
This is explained by the circularizatio n by t ides raised on 
the star by the planet (Marcy et al. 1997 ). One excep- 



tion to this rule, HD83443b (e = 0.079 ± 0.033), can be 
attributed to the presen ce of another eccentric planet in 
the system (Mayor et al. 2001 ) . As we shall see, the plan- 
ets inside ~0.1 AU also have very specific properties due 
to the closeness to their star and the intense radiation 
they receive. For this reason, following astronomical con- 
ventions, we choose to name them after the first object of 
this class to have been discovered: "51 Pegb-like" planets, 
or in short "Pegasi planets" . 

Such planets provide an unprecendented opportunity 
to study how intense stellar irradiation affects the evo- 
lution and atmospheric circulation of a giant planet. 
Roughly 1% of stars surveyed so far bear Pegasi planets in 
orbit, suggesting that they are not a rare phenomenon. 
Their proximity to their stars increases the likelihood that 
they will transit their stars as viewed from Earth, allow- 
ing a precise determination of their radii. (The probability 
varies inversely with the planet's orbital radius, reaching 
~10% for a planet at 0.05 AU around a solar-type star.) 
One planet, HD209458b, has already been observed to 



transit its star every 3.524 days (Charbonncau et al. 2000 



Henry et al. j2000[ ). The object's mass is 0.69 ± 0.05 Mj , 
where Mj= 1.89 x 10 27 kg is the mass of Jupiter. Hubble 
Space Telescope measurements of the transit (Brown et al. 
|200lD imply that the planet's radius is 96300 ± 4000 km. 
An analysis of the lightcurve combined with atmospheric 
models shows that this should correspond to a radius of 
94430 km at the 1 bar level (Hubbard et al. |200l| ). This last 
estimate corresponds to 1.349 Rj, where Rj = 70,000km 
is a characteristic radius of Jupiter. This large radius, in 
fair a greement with theoretical predictions (Guillot et al. 



1996 ), shows unambiguously that HD209458b is a gas gi- 



ant. 

We expect that the evolution of Pegasi 
planets depends more on the stellar irradiation than 
is the case with Jupiter. HD209458b and other 
Pegasi planets differ qualitatively from Jupiter be- 
cause the globally-averaged stellar flux they absorb is 
~ 10 8 ergcm~ 2 (10 5 Wm~ 2 ), which is ~ 10 4 times greater 
than the predicted intrinsic flux of ~ 10 4 ergcm~ 2 . (In 
contrast, Jupiter's absorbed and intrinsic fluxes are the 
same within a factor of two.) Several evolution calcula- 
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tions of Pegasi planets have been published (Burrows et 
al. 2000a, Bodenheimer et al. |200l| ), but these papers 
disagree about whether HD209458b's radius can be 
explained, and so far there has been no general discussion 
of how the irradiation affects the evolution. Our aim is to 
help fill this gap. 

Here, we quantify how atmospheric processes affect the 
evolution of Pegasi planets such as HD209458b. First (Sec- 
tion 2), we show that the evolution is sensitive to the as- 
sumed atmospheric temperatures. This sensitivity has not 
previously been documented, and quantifying it is impor- 
tant because the temperature profiles appropriate for spe- 
cific planets remain uncertain (e.g., no atmospheric radia- 
tive transfer calculation for HD209458b yet exists). Our 
works suggests that the discrepancy between the predic- 
tions of Burrows et al. ( [2000a ) and Bodenheimer et al. 



(2001) can be largely explained by their different assump- 
tions about atmospheric temperature. 

Second, the effect of atmospheric dynamics on the evo- 
lution has to date been neglected. For example, current 
models assume the day-night temperature difference is 
zero, despite the fact that substantial day-night temper- 
ature variations are likely. In Section 3.1 we demonstrate 
how the evolution is modified when a day-night tempera- 
ture difference is included. Furthermore, the intense stellar 
irradiation will lead to production of atmospheric kinetic 
energy, and transport of this energy into the interior could 
provide a substantial energy flux that would counteract 
the loss of energy that causes planetary contraction. In 
Section 3.2 we investigate this effect. 

The research has major implications for HD209458b. 
Early calculations implied that Pegasi planets contract 
slowly en ough to explain HD20 9458b' s large radius (Guil- 



lot et al. 1996| , Burrows et al. 2000a ). But recent calcu- 



lations of irradiated atmospheres suggest that the actual 
deep atmosphere is colder than assumed (Goukenleuque et 
al. 2000). When such realistic temperatures are adopted 
(our Section 2), the planet contracts too fast and t he ra - 
dius is ~ 0.2-0.3 Rj too small. Bodenheimer et al. ( [2001 ) 
argued that tidal heating from circularization of the orbit 
would slow the contraction, leading to a larger radius, but 
this is a transient process that would end ~ 10 8 years af- 
ter the planet's formation. Instead we argue that kinetic 
energy produced in the atmosphere is transported into the 
interior and dissipated (Section 3.2). We show that plau- 
sible downward energy fluxes can slow or even halt the 
planet's contraction, allowing HD209458b's radius to be 
explained. 

In a joint paper (Showman & Guillot [2001 , Paper II) 
we consider the atmospheric dynamics of these planets, 
with emphasis on how the atmospheres respond to stellar 
heating and gravitational tidal interactions, and on the 
observable consequences. 



2. Sensitivity of Evolution to Atmospheric 
Temperature 

The upper boundary condition of evolution models con- 
sists of a relationship between the effective temperature 
and some deeper temperature (say that at 10 bars) to 
which the model's interior temperature profile is attached. 
Here we show that the evolution is sensitive to the as- 
sumed relationship (i.e., to the assumed atmospheric tem- 
perature structure). 

Before we begin, we provide some definitions. We de- 
fine the effective temperature of the irradiated planet as 

AnR 2 (jT^ = L 9 +L int , (1) 

where R is the planet's radius, T e fj is its effective temper- 
ature, cr is the Stefan-Boltzmann constant, L 9 is the part 
of the stellar luminosity absorbed by the planet and L; nt is 
the intrinsic luminosity of the planet due to its cooling and 
contraction (and possibly other processes such as radioac- 
tivity or thermonuclear reactions in the case of massive 
objects). 

The temperature corresponding to the intrinsic plane- 
tary flux, called the "intrinsic" temperature Ti nt , is defined 

by 

^R 2 aT^ t = L int . (2) 

Similarly, we define T from the absorbed stellar lumi- 
nosity L 9 . T 9 is the effective temperature towards which 
the plane t tend s as it cools and L lni diminishes (see e.g. 
Hubbard 1977 ). It is a function of the Bond albedo (i.e. 
the ratio of the luminosity directly reflected to the total 
luminosity intercepted by the planet). T a can be viewed 
as the effective temperature reached by the planet after it 
has lost its internal heat, and is hence sometimes noted 



T eq (e.g. Guillot et al. |1996| ; Saumon et al. |1996[ ). 



2.1. Atmospheric boundary conditions 

We consider two evolution models of HD209458b based 
on the parameters listed in table [j]; the two models differ 
only in their prescription for the atmospheric boundary 
condition. 

Our first evolution sequence, dubbed the "hot" case, 
uses the standard boundary condition from Guillot et al. 
(1996) and Burrows et al. (2001a). These papers adopted 
an atmospheric structure of an isolated object with the 
expected effective temperature, which provides a fair fit 
to the evolution of Jupiter. The surface boundary condi- 
tion consists of a relationship between the temperature at 
the lObar pressure level Ti so i atc d, the effective tempera- 
ture T e g and the gravity g of an isolated planet/brown 



dwarf derived by several authors (see Mar ley et al. 1996 
Burrows et al 



1997) 



T(P = 10 bar) = T iso \ at ed{T e s , g) 



(3) 



This approximation is exact in the limit when the stellar 
flux is entirely absorbed at the 10 bar level, or in a deep 
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Parameter 


Value 


References /remarks 


Evolution model 


CEPAM 


Guillot & Morel (1995) 


Mass 


M 


= 0.69 M Jup 


(M Jup = 1.89 x W w g) 


Absorbed stellar heat 


T 9 


= 1400 K 




Radius 


R- 


= 1.35 R.j 


(Rj= 7 x 10 9 cm) 



EOS "interpolated" 

Helium mass mixing ratio Y = 0.30 

Rosseland opacities 

Rotation 

Core mass 

Atmospheric boundary 



Saumon, Chabrier & Van Horn (1995) 

Higher than solar in order to mimic a solar abundance of 

heavy elements 

Alexander & Ferguson (1994) inch interstellar grains 
Neglected 
Not considered 



From Marley et al_( [l996D ; Burrows et al. ( |l997l ) 
See Eqs. <Jsb and " 
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Fig. 1. Surface boundary condition (temperature at the 
10 bar level) that has been used in several published evo- 
lution models, and which we dub the "hot" case, as a 
function of effective temperature for three different gravi- 
ties: log(g)=2.5, 3 and 3.5. (g is in units of cms -2 ). (See 
Saumon et al. (1996) for a discussion.) 



adiabatic (convective) region, as is the case for Jupiter 
(see Hubbard |1977 for a detailed discussion of the effect 
of insolation on Jupiter's evolution). Figure [l] shows the 
variation of Isolated with effective temperature and grav- 
ity, for values of interest in the case of Pegasi planets. 

Unfortunately, the approximation becomes incorrect 
in the case of strongly irradiated planets because of the 
growth of a thick external radiative zone. Another bound- 
ary condition has therefore to be sought: either part 
of the stellar flux is able to penetrate to deeper levels 
(Po > 10 bar) and lead to a boundary condition defined 
by T(P ) > Isolated, or most of the stellar flux is absorbed 
at Po < 10 bar, yielding T(P ) < Ti so i ate d. (This is due to 
the fact that in the radiative zone dT/dP oc F, where F 
is the flux t o be transported). It will be shown hereafter 
(see section 3J3) that Eq. ||| is effectively an upper limit 
to the boundary temperature because continuum opacity 



sources only effectively limit the penetration of the stellar 
photons. Indeed, more detailed models of the atmospheres 
of Pegasi planets have shown that most of the starlight is 
absorbed at pressures less than 10 bar, and that Eq. || 
overestimates the atmospheric temp eratu r es by as much 
as 300 to 10 00 K ( Seager & Sasselov |l998l £000l Gouken- 



leuque et al. 2000 ; Barman et al. 2001 ). 

Because these atmospheric models do not presently 
span the effective temperature and gravity range that is 
needed, and more importantly because they assume un- 
realistic intrinsic temperatures, we chose to construct an 
arbitrary boundary condition based on the results of the 
isolated case. For a given T e g, the isolated case provides 
an upper bound to the "surface" temperature and by ex- 
tension to the temperatures in the planetary interior. In 
order to have an approximate lower bound that agrees 
with atmospheric models of irradiated giant planets, we 
assume (i) a lower value of Po = 3 bar, and (ii) that the 
temperature at that level is given by: 



T(P = 3 bar) = T isolatod (T eff , g) - 1000 K. 



(4) 



Evolution calculations made with this boundary condition 
are dubbed the "cold" case. 

Note that we found a posteriori that the choice of Po 
is almost unconsequential for the evolution calculations. 
This is because the external radiative region quickly be- 
comes almost isothermal (see fig. || hereafter). However, 
the consequences of the cooler temperatures are profound, 
and as we shall see lead to a much faster evolution. 



In this context, Bodcnhcimcr et al. (2000; 2001) as- 
sume that the temperature at optical depth 2/3 (cor- 
responding in their model to a pressure of the order 
~ 1 mbar) is equal to the effective temperature T c g , an 
approximation that leads to an underestimation of the 
actual atmospheric temperatures. As a consequence, their 
1 bar temperatures are of the order of ~ 1400 K, i.e. even 
lower than what is implied by Eq. (0). This would imply 
an extremely inefficient penetration of the stellar flux in 
the planetary atmosphere, in disagreement with detailed 
models of these atmospheres. We therefore prefer to use 
Eq. (0) as our "cold" boundary condition. 
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2.2. Evolution models of HD209458b 

The evolution of HD209458b is calculated as described in 



Guillot ct al. (1996), using the parameters given in Ta- 
ble 0. Because of the high stellar insolation, the contrac- 
tion and cooling of the planet from a high entropy initial 
state is only possible through t he bu ild-up and growth of 
a radiative zone (Guillot et al. 



19961) 




1.0 8.5 9.0 
Log(Time) [yr] 



10.0 



Fig. 2. Evolution of HD209458b using the "hot" atmo- 
spheric boundary condition (Eq. (||)). The evolution of 
the central pressure with time is shown as the bottom 
thick line. The planet is convective except for an upper 
radiative zone indicated by a hashed area. Isotherms from 
4000 to 20 000 K are indicated. The isotherms not labelled 
correspond to 3500, 30 000 and 40 000 K. The dashed line 
indicates the time necessary to contract the planet to a 
radius of 1.35 Rj. [1 bar = 10 6 dyn cm . 



The evolution of the interior of HD209458b for the two 
cases is shown in figures |2| and ||. After a rapid contrac- 
tion during which both the central pressure and temper- 
ature increase, the onset of degeneracy leads to a cooling 
of the interior as the planet continues to contract. The 
cooling of the interior proceeds despite the fact that the 
atmospheric temperatures remain nearly constant thanks 
to the growth of a radiative region in the planet's upper 
layers, as indicated by the dashed area. 

In the case of the "hot" atmospheric boundary condi- 
tion (Fig. the measured radius (1.35 Rj) is attained af- 
ter 5.37 Ga (5.37x 10 9 years), which is compatible with th e 
age of the GO star HD209458 (see Burrows et al. |2000a|) . 
The radiative zone then extends to about 730 bar, and the 



intrinsic luminosity is 7.7 x 10 24 ergs 1 (2.3 times that of 
Jupiter), which corresponds to an intrinsic temperature of 
105K. 

The "cold" atmospheric boundary condition (Fig. |J) 
yields a much faster evolution: the planet then shrinks to 
1.35 Rj in only 0.18Ga (see also fig. ^ hereafter). This is 
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Fig. 3. Same as Fig. || but with the "cold" boundary 
condition given by Eq. (||). Unlabeled isotherms are for 
T = 2500, 30 000 and 40 000 K. 



incompatible with the age derived for HD209458b. The 
radiative/convective boundary in this model (at 0.18 Ga) 
is at 160 bar, due to the higher intrinsic luminosity equal 
to 1.9x 10 26 ergs _1 , equivalent to an intrinsic temperature 
of 234 K. 

The fact that the evolution is faster in the "cold" case 
may seem counterintuitive. It occurs because the intrinsic 
luminosity is proportional to the temperature gradient, 
not to the temperature itself. As shown by Fig. [i] hereafter, 
the temperature variation in the radiative region is more 
pronounced in the "cold" case than in the "hot" case. 
Basically, this is because the temperatures at deep levels 
are fixed by the condition on the radius, but that the 
surface temperatures are very different in the "cold" and 
"hot" cases. 

Bodenheimer et al. (2001) obtain radii that slightly 
exceed those of our "cold" case, despite their lower at- 
mospheric temperatures. However, this is probably due 
to their lower assumed value for the helium abundance 
Y = 0.24, whereas we chose Y — 0.30, a value represen- 
tative of conditions in the solar nebula and that accounts 
for a solar proportion of heavy elements. In both cases, 
young ages are required to reproduce the measured plan- 
etary radius. 

We therefore feel that because the "cold" atmospheric 
boundary condition is preferable to the "hot" boundary 
condition, there is a problem in explaining HD209458b's 
radius. An absolute proof of this statement would require 
calculations of many different models using different as- 
sumptions, which we will not attempt in this paper. The 
conclusion should be relatively secure however, because 
several factors point towards a reduction of the planet's ra- 
dius compared to what we have calculated: (i) The atmo- 
spheric temperatures could be even lower than in the cold 
case; (ii) The opacities used include the presence of abun- 
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dant grains in the atmosphere and does not account for 
their gravitational settling; (iii) Our choice of the equation 
of state tends to yield larger radii than would be the case 
using an equation of state that consistently models the 
molecular/metallic transition (see Saumon et al. |1995| ); 
(iv) The presence of a central core will tend to greatly 



reduce the planet's radius (Bodcnhcimcr et al. 2001). 

An additional source of energy then appears to be re- 
quired. We note that the presence of a hydrogen/helium 
phase separati on, lik e in Jup iter a nd Saturn (e.g. Steven- 
son & Salpeter [l97^ ; Guillot |l999|) , is not a valid alterna- 
tive because of the high interior temperatures involved in 
the case of Pegasi planets (e.g. 20,000K at IMbar). 

An important aspect of the hot and cold models is 
shown in Fig. ^: Apart from the outer radiative layers, 
the two models possess a very similar interior structure at 
the times (5.37 and 0.18 Ga for the hot and cold cases, 
respectively) when they match HD209458b's observed ra- 
dius. This is easily understood by the fact that the radia- 
tive layer encompasses only a small fraction of the radius. 
Most of the contribution to the planetary radius is due to 
the convective interior. Fixing the radius is, for a given 
equation of state and composition, almost equivalent to 
fixing the temperature-pressure profile in the deep inte- 
rior. HD209458b can be thought of as a relatively well- 
constrained convective core underlying a radiative enve- 
lope of uncertain mass and temperature. 
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3.1. Implications of atmospheric day/night temperature 
variations 

We show in Paper II that the atmosphere should be sig- 
nificantly hotter on the dayside than the night side. Here 
we examine the consequences for the evolution models. 
To do so, we use a toy model that, while simple, eluci- 
dates the important physics. A full two-dimensional model 
that would allow us to calculate the evolution of Pegasi 
planets including latitudinal or longitudinal temperature 
variations will be left for future work. 

Let us assume that the planet can be divided in two 
hemispheres (night and day) with two different effective 
temperatures such that Tnight < ?day- When the absorbed 
stellar energy is fully redistributed by advection, Tday = 
Tnight = T 9 . In all cases, energy conservation implies that 
2T£ = T* ight + T^ ay . Therefore, T night = implies T da y = 
2 i/4 Ta . Ta = 1400K and T day = 1500 K yields T night = 
1272 K. 



Stellar flux 




Fig. 5. Schematics of the day-night toy evolution model. 
The slow mixing of the interior leads to a non-radial heat 
flux Lday-night f rom the day side to the night side. As a 
consequence, the intrinsic luminosity on the day side is 
smaller and that on the night side becomes larger. 
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Fig. 4. Temperature profiles for the 5.37Ga-old "hot" 
model (thin black line) and the 0.18Ga-old "cold" model 
(thick grey line), i.e. when they match HD209458b's 
measured radius. The diamonds indicate the radia- 
tive/convective boundary. 



3. Evolution of Pegasi planets: "non-standard" 
models 



Let us do the following gedanken experiment, as illus- 
trated by Fig. ^: we suppose that the two hemispheres can- 
not exchange energy, and let them evolve from the same 
initial state. After a given time interval At, the central 
entropy on the day side will have decreased by a smaller 
amount than on the night side. This is due to the fact that 
a higher atmospheric temperature is equivalent to a higher 
stellar flux, and leads to a slower evolution (see Hubbard 



1977 ; Guillot et al. 1995 ). In consequence, the night side 



will have become internally colder, have a smaller radius 
and a larger central pressure than the day side. 

The pressure differences caused by the differential cool- 
ing ensures an efficient mixing between the two hemi- 
spheres on a time scale of decades or less, i.e. much shorter 
than the evolution time scale. 
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We therefore include the effect of atmospheric temper- 
ature variations on the evolution in the following way: We 
calculate two evolution tracks of a planet with uniform 
temperatures T„i g ht and Td ay , respectively. Using these 
evolution tracks and starting from an initial condition for 
which the two models have the same central entropy, we 
calculate the entropies of the two sides after a time inter- 
val St. We then decrease the entropy of the day side and 
increase that of the night side so that both are equal to 
(knight + 5 l c j a y)/2. The process is repeated for each time 
step. The cooling of the night side is therefore slowed by 
the mixing of material of slightly higher entropy from the 
day side, while the opposite is true for the day side. 



given mass, composition and stellar insolation, up to 5% 
smaller than calculated in the uniform case. 
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Fig. 6. Radius of HD209458b (in units of the radius of 
Jupiter) versus time under different assumptions. The 
plain lines correspond to the "hot" and "cold" evolution 
cases shown in figures ^] (upper curve) and |3| (lower curve), 
respectively. The long dashed line is obtained in the "cold" 
case, when assuming that the radiative equilibrium effec- 
tive temperature is 1500 K on the day side and 1272 K on 
the night side. The short dashed line is obtained when 
these effective temperatures are 1664 K and OK, respec- 
tively (~ no advection). The box indicates inferred radii 
and ages of HD209458b. 



The resulting evolution tracks are shown in Fig. ^. 
Not suprisingly, the cooling of an irradiated planet with 
inefficient heat redistribution in the atmosphere is faster 
than if the stellar heat is efficiently advected to the night 
side. This is mainly due to the fact that, with increasing 
AXday-night, T n i gn t decreases much more rapidly than Tday 
increases, yielding a much faster cooling of the night side. 
However, for the temperature variation of 228 K shown in 
Fig. U (long dashed line), the effect is limited to a variation 
of ~ 0.5% of the radius after 1 Ga of evolution or more. 
The effect is of course more pronounced if no thermal en- 
ergy advection occurs in the atmosphere (T n ; g ht = 0). In 
that case (short dashed line), the minimal radius is, for a 



Fig. 7. Intrinsic planetary fluxes obtained as a function 
of time in the "cold" case with a well mixed atmosphere 
(plain line) . When assuming that stellar irradiation is im- 
perfectly redistributed over the planet's atmosphere, the 
flux on the night side (Tnight =1272 K) becomes larger 
than that on the day side (Td ay =1500 K). 



Figure shows that a non-uniform atmosphere 
has a substantial effect on the planet's intrinsic flux 
(F = L/AttR 2 ). After 4.5 Ga, the flux of our planet 
with T night =1272 K and T day =1500 K is 6310 and 
2950 erg cm -2 s -1 on the night and day sides, respec- 
tively which can be compared with the measured in- 
trinsic fluxes of Jupiter, 5440 erg cm -2 s -1 and Saturn, 



2010 erg cm s (Pearl & Conrath 1991). This process 



is analogous to that proposed by Ingersoll & Porco (1978) 



to explain the uniform temperatures of Jupiter. Stellar in- 
solation tends to suppress the planet's intrinsic heat flux, 
and so the planetary heat preferentially escapes in regions 
where the insolation is minimal. 



3.2. Evolution with internal energy dissipation 

Current models predict that several-Ga-old Pegasi 
planets have intrinsic heat fluxes of ~ 10 4 erg s -1 cm -2 , 
which is about 10 4 times less than the total luminos- 
ity of ~ 10 s erg s -1 cm -2 resulting from thermal bal- 
ance with the stellar insolation. A fraction r\ of the 
total luminosity will be converted into kinetic energy 
by the atmospheric pressure gradients. On Earth, the 
globally-averaged flux transported by the atmosphere is 
about 2.4 x 10 5 ergs -1 cm -2 (240 Wm -2 ), while about 
2000 ergs -1 cm -2 is converted into large-scale atmo- 



spheric kinetic energy (Peixoto and Oort 1992, p. 385), 
leading to a value r\ — 0.01. This energy production can 
be viewed as the work done by an atmospheric heat engine 
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with an efficiency of 1%. Preliminary simulations that we 
have conducted indicate that a similar ratio is relevant for 
Pegasi planets (Paper II). If so, the implied kinetic energy 
generation is 10 2 times the intrinsic heat flux computed 
by current models. Inclusion of this energy could then lead 
to a first-order alteration in the behavior predicted in evo- 
lution models. 

In steady state, the kinetic energy that is produced 
must be dissipated. On Earth, this dissipation mostly re- 
sults from friction with the surface (Peixoto and Oort 
1992, p. 385). For Pegasi planets, Kelvin-Hclmholtz in- 



stabilities and breaking of gravity and planetary waves 
are more relevant. The key question is whether the en- 
ergy is dissipated in the "weather" layer, where starlight 
is absorbed and radiation to space occurs, or in the deeper 
atmosphere. In the former case, the dissipation will pro- 
vide only an order- 1% perturbation to the vertical profile 
of absorbed starlight and radiated thermal energy. In the 
latter case, it comprises a hundred-fold alteration in the 
interior energy budget. We therefore need to know (i) to 
what pressures can the energy be transported, and (ii) 
how deep must it be transported to cause a major effect 
on the interior? 

As discussed in Paper II, mechanisms for transporting 
kinetic energy into the interior include Kelvin-Hclmholtz 
instabilities, direct vertical advection, and waves. The 
dynamical coupling between atmospheric layers suggests 
that winds should develop throughout the radiative re- 
gion even though the radiative cooling and heating oc- 
curs predominantly at pressures less than a few bars. The 
boundary between the radiative region and the convective 
interior (at 100-1000 bars depending on the model) is a 
likely location for dissipation, because Kelvin-Hclmholtz 
instabilities and breaking of downward propagating waves 
can both happen there. Furthermore, application of the 
Taylor-Proudman theorem to the convective interior sug- 
gests that winds should develop throughout the convective 
interior even if the forcing occurs only near the top of the 
convective region. This increases the possibility of dissipa- 
tion in the interior. 

With the inclusion of an internal dissipative source, 
the energy equation becomes 



dL_ 

dm 



T dS 
dt' 



(5) 



where m is the mass inside any given level, and e(m) is 
the energy dissipated per unit time per unit mass at that 
level. 

The evolution of Pegasi planets including energy dis- 
sipation has been studied by Bodenheimer et al. ( |2001 ) 
in the context of the tidal circularization of the orbit of 
the planet. These authors focused on simulations where e 
was constant with m (although they also performed some 
simulations with spatially- varying dissipation). The ma- 
jor difficulty is, as noted by the authors, the fact that 
the present eccentricities of extrasolar planets within 0.1 



AU of their star are small and that the detected Pegasi 
planets generally do not possess close massive planetary 
companions which would impose on them a forced eccen- 
tricity. 

Instead, we have argued that kinetic energy, gener- 
ated from a portion of the absorbed stellar flux, is trans- 
ported to the interior where it can be dissipated. Although 
the depth of such dissipation is unknown, the majority 
could be deposited within the radiative zone rather than 
throughout the interior. Due to the rapid rise of the Rosse- 
land opacity with pressure and temperature, the effect of 
heating anywhere within the convective core is essentially 
equivalent to the case where it occurs entirely at the cen- 
ter (a result shown by Bodenheimer et al.). The question 
is whether even shallower heating — say that occurring 
at tens to hundreds of bars, where atmospheric kinetic- 
energy deposition is likely — can affect the evolution. 
Therefore, we here explore the influence of the dissipa- 
tion's depth dependence and magnitude E = J edra. 

An ad hoc, but reasonable, assumption is that a frac- 
tion of up to 1% of the absorbed stellar flux is dis- 
sipated inside the planet. Quantitatively, we use E = 
2.4 x 10 27 ergs _1 . Relatively small values of E can affect 
the evolution, provided they are comparable or larger than 
the luminosity obtained without dissipation L (note that 
L ~ 10 24 — 10 25 ergs _1 ) and affect the radiative gradient 
on a sufficiently extended region of the interior. 

A first calculation assumes that energy is dissipated 
entirely at the center of the planet. In that case, as shown 
in Fig. ^| (uppermost solid line), an equilibrium with the 
star is reached after only ~ 100 Ma, at which point the 
planet's radius is 1.87 Rj and its structure remains un- 
changed with time (as long as the star is also in equi- 
librium). This is very similar to the results obtained by 
Bodenheimer et al. (2001). Also in agreement with their 
results, we find that a calculation with the same E, but 
with the dissipation evenly distributed throughout the in- 
terior (i.e., uniform e) yields a curve similar to the upper 
solid line in Fig. |§|. 

In order to estimate the consequences of energy dissi- 
pation occurring closer to the planet's surface, we use the 
following arbitrary functional form: 



(6) 



where m is the adimensional mass (0 at planet's center and 
1 at its surface), mo is the mass fraction of the external 
regions over which most of the dissipation occurs, and eo 
is chosen such that J edm = E = 2.4 x 10 27 ergs _1 . We 
will use values of mo equal to 5 x 10~ 6 , 10~ 5 , and 2 x 10~ 5 , 
which correspond to locations where the pressure is 5, 11, 
and 21 bars, respectively. 

A choice of mo = 10~ 5 (which implies dissipative 
heating distributed dominantly from the top boundary 
to 15 bar but with a tail of heating reaching ~ 100 bar) 
yields a radius which is, after a few billion years, about 
10% larger than in the case with no dissipation (third solid 
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Fig. 8. Evolution tracks obtained in the "cold" case, show- 
ing the influence of dissipation. The bottom grey line cor- 
responds to the case with no dissipation. The other solid 
lines have been calculated including the dissipation of 1% 
of the absorbed stellar flux (2.4 x 10 27 ergs _1 ), at various 
depths: from bottom to top, dissipation was supposed to 
occur in various mass shells: mo = 10~ 5 , 2 x 10 -5 , or at 
the center of planet, respectively. The dashed line corre- 
sponds to dissipation of 10% of the absorbed stellar flux 
(2.4 x 10 28 ergs _1 ) with m = 5x 1CT 6 . The dotted evolu- 
tion track is for dissipation of E = 1.8 x 10 26 erg s _1 at the 
planet's center. In the first two cases, energy dissipation 
occurs mostly from the upper boundary to the nearly- 
isothermal region. The m = 5 x 10~ 6 , 10~ 5 and 2 x 10~ 5 
values of adimensional mass correspond to pressures of 5, 
11, and 21 bars, respectively. 



line from the top in Fig. ||). This is, in the "cold" case, in- 
sufficient to reproduce the observed radius of HD209458b. 
A slightly higher value of too = 2 x 10~ 5 yields an evo- 
lution track which is in agreement with the measured ra- 
dius, as shown in Fig. ^ (second solid line from the top). 
In that case, the value of 1 — rh = too corresponds, for 
the model with 1.35 Rj to a pressure level P = 28 bar and 
T = 2800 K. However, because of the form of Eq. (0), dis- 
sipation becomes negligible only around 1 — rh w 10mo, 
i.e. P = 130 bar and T = 3380 K. 

Two other evolution tracks have been calculated 
specifically to illustrate how HD209458b's radius can be 
reproduced with different values of E and mo- In the case 
of dissipation at the center, we were able to match an equi- 
librium radius of 1.35 Rj with E = 1.8 X 10 26 ergs _1 (dot- 
ted line in Fig. |8|), which is only 0.08% of the global-mean 
absorbed stellar flux. In the case of dissipation limited to 
a shallow layer (too = 5 x 10 -6 , corresponding to a pres- 
sure of 5 bars), we found that a relatively high value of 
E corresponding to 10% of the absorbed stellar flux was 
necessary for the planet to contract to its present radius 
in about ~ 5 Ga (dashed line in Fig. ||). 



Fig. 9. Temperature-pressure profiles for models of 
HD209458b (with R = 1.35 Rj) calculated with the "cold" 
atmospheric boundary condition, and assuming various 
dissipation profiles: E = 1.8 x 10 26 ergs _1 at the planet's 
center (solid line), E — 2.4 x 10 27 ergs _1 (1% of the 
globally-averaged absorbed stellar heat flux) and too = 
2 x 10~ 5 (dashed line), E — 2.4 x 10 28 ergs _1 and rh = 
5 x 10~ 6 (dot-dashed line). The diamonds indicate the ra- 
diative/convective boundary. In the case of the dot-dashed 
models, the large dissipation in the external layers is re- 
sponsible for a second convective zone pressures lower than 
25 bars, as indicated by a triangle. 



The temperature pressure profiles of three models cal- 
culated with the "cold" atmospheric boundary conditions 
but with different values of the dissipation factor, and such 
that their total radius is 1.35 Rj , are compared in Fig. ^. 
The temperature profile of the model with dissipation at 
the center (solid line) is essentially indistinguishable from 
our reference "cold" model which included no dissipation 
but had a high intrinsic heat flux due to its young age. 
As dissipation is increased but at the same time limited 
to shallower outer layers, the temperature profile becomes 
more similar to the "hot" case shown in Fig. |]. In the case 
of the highest dissipation considered here but with a small 
too = 5 x 10~ 6 , an external (detached) convective region 
can form. Note that in this case Tint = 790 K; atmospheric 
models calculated with these high intrinsic effective tem- 
peratures are also found to possess a deep convective zone 
(Barman et al. |200l|). 



3.3. Penetration of the stellar 
atmosphere 



lux into the deep 



We have shown that a kinetic energy flux corresponding 
to a small fraction of the stellar flux can, if dissipated 
deep enough, significantly affect the planet's evolution. 
The result would be exactly the same if the stellar flux 
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was radiatively transported to these deep levels. We argue 
that stellar heat cannot be deposited so deep, however. 

As shown by figures [| and ||, the temperature profile of 
HD209458b must cross the point defined by P w f kbar, 
Ti « 4000 K, assuming that the planet is of solar composi- 
tion (the addition of heavy elements tends to increase the 
temperature required at a given pressure, but doesn't oth- 
erwise alter the conclusions that follow). This point and 
the external boundary condition then define the intrinsic 
luminosity required to reproduce the measured radius, i.e. 



H 2 -H 2 collision-induced absorption (see e.g. Lenzuni et al. 
Guillot ct al. 



MnaGM < T 3 > / dT 
3 < k > \ dp 



(7) 



rati 



where a is the Stefan-Boltzmann constant, < T > and 
< k > are a mean temperature and opacity, respectively, 
and (dT/dp) ra( j is the mean temperature lapse rate in the 
radiative zone. 

In the "hot" case, the value of Li nt thus derived is small 
because the difference in temperature between the bottom 
of the radiative zone and the external boundary is small. 
This implies that HD209458b needs to be relatively old to 
have such a low intrinsic luminosity. In the more realistic 
"cold" case, the planet either has to be uncomfortably 
young, or some additional heat has to be transported to 
these levels. This requires e > 0, but with the additional 
requirement that the temperature should be brought close 
to ~ 3200 K at a pressure P < p. 

By definition of the optical depth t\, the proportion 
of stellar flux still remaining at a given level is equal to 



where 



K\pdz, 



(8) 



and n\ is the monochromatic opacity, and z is the altitude 
in the atmosphere. Equation (|^) is strictly valid only in a 
one-dimensional approximation, but this is sufficient here. 
It is useful to approximate the integral by assuming a 
constant opacity and density in a slice of atmosphere of 
height equal to the pressure scale height. In that case, the 
pressure at which 99% of the stellar flux has been absorbed 
can be estimated by: 



P £ 5- 
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(9) 



where K m i n corresponds to the minimum value of the opac- 
ity K\. 

We estimate from the previous section that a pene- 
tration of 1% of the stellar flux to P ~ 100 bar in the 
"cold" case allows the radius of HD209458b to be ex- 
plained without any other energy dissipation. 99% of the 
flux of a 6000 K black body is emitted between 0.22 and 
4.9 /im. The measured radius and mass implies that g w 
10 3 cms -2 , therefore requiring K m ; n < 5 x 10~ 5 cm 2 g _1 . 

This opacity is approximately the minimum expected 
for a pure hydrogen-helium mixture at P ~ 1 bar and T ~ 
1500 K, at A ~ 1 /im due to Rayleigh scattering by H 2 and 



1991 



1994). At temperatures above about 



2000 K, two very important sources of continuous opacity 
arise, led by the increasing number of free electrons: the 
free-free absoption of H^ and the bound-free absorption 
of the H~ ion. However the number of free electrons in a 
zero-metallicity gas remains low even at 3000 K, and the 
low opacity minimum persists to temperatures exceeding 
3000 K, and pressures exceeding 10 bar. In this case a deep 
absorption of the stellar flux would then be likely. 

However, in a mixture of solar-like composition, a large 
fraction of the electrons can be provided by alkali met- 
als. Using electrons number densities obtained from Ku- 
rucz ( 1970 ) and Lodders (personal communication, 2001), 
we estimate the minimum continuous opacity to climb to 
3 x 10- 3 cm 2 g- 1 at 2500 K and to 0.1 cm 2 g^ 1 at 3500 K, 
mostly due to H~ absorption (John 1988 ). This alone pre- 
vents any relevant fraction of the stellar flux to reach levels 
at which the temperature is larger than 2500 K. 

Furthermore, a number of other opacity sources are ex- 
pected to occur and even dominate the spectrum. Likely 
candidates are K and Na which are now known to con- 
tribute significantly to the atmospheric absorption of 
brown dwarfs with similar temperatures, at visible wave- 
lengths (Burrows et al. [20001 ). Similarly, TiO is expected 
to provide an even larger absorption at short wavelength 
where it appears in the deeper atmosphere. For example, 



Barman et al. ( |2001 ) find that ri^m = 10 is attained 



at pressures smaller than ~ 6 bar in the cloud-free atmo- 
sphere of Pegasi planets. Finally, clouds, if present, would 
cause an absorption of the stellar flux at even lower pres- 
sures. 

It hence appears that only a zero-metallicity atmo- 
sphere would have a low-enough opacity to allow the stel- 
lar flux to penetrate to P ~ 100 bar. This is an unlikely 
possibility, the metallicity of HD209458 being close to so- 
lar (Mazeh et al. 20001 ). One possibility remains however: 
that alkali metals and strong absorbers such as TiO are 
buried deep due to condensation effects on the night side 
(see Paper II), so that the atmosphere on the day side 
would be almost metal-free. It is not clear even in this case 
that the measured radius could be explained, because the 
lower overall opacities would increase the rate of cooling 
and hence contraction of the planet. 

In all the cases considered here, it seems very difficult 
for the incoming stellar flux to penetrate down to levels 
where the temperature is large (more than ~ 2500 K). 
In order to reproduce HD209458b's large radius, a tem- 
perature ~ 4000 K at a pressure P ~ 1 kbar must be 
attained. Energy dissipation due to a transfer of kinetic 
energy hence appears as the most likely missing energy 
source. 
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4. Conclusions 

We have shown that the evolution of Pegasi planets is 
mainly driven by processes occuring in their atmosphere 
and is consequently complex. The measurement of the ra- 
dius of one of these objects, HD209458b, has allowed us 
to probe some of these mechanisms in detail. 

We demonstrated that radiative-equilibrium atmo- 
spheric models predicting temperatures above ~ 2500 K 
at pressures P ~ 10 bar are unlikely given the rapid rise of 
the absorption with increasing temperature. Cooler tem- 
peratures are to be expected in the atmosphere and with- 
out other means than radiation to transport the incoming 
heat flux, HD209458b's large radius cannot be reproduced 
unless the planet is much younger than is revealed by ob- 
servations of its parent star. 

We showed the atmospheric temperature variations to 
have a small effect on the planetary cooling, if limited to 
a few 100's K. The temperature variations lead to faster 
cooling of the planet compared to standard models, which 
assume the stellar heat to be evenly distributed onto the 
planet's atmosphere. This accentuates the problem of re- 
producing HD209458b's radius. 

Energy dissipation is however a very promising can- 
didate to explain HD209458b's missing heating source. 
Lubow et al. (1997) have shown that tidal synchroniza- 
tion of Pegasi planets could give rise to a large heat flux. 
But this mechanism is limited to the early evolution of 
the planet and should rapidly become negligible. Boden- 
heimer et al. (2001) argued that internal heating could be 
provided by tidal circularization of an eccentric orbit. This 
is similarly unlikely to occur in most Pegasi planets in the 
absence of a detected close, massive companion capable of 
exciting their eccentricity. The mechanism that we invoke 
is simply a downward transport of kinetic energy gener- 
ated by the intense atmospheric heat engine. We showed 
that only ^0.08% of the stellar flux has to be transported 
to the interior regions to explain the radius of HD209458b. 
This fraction rises to 1% if heat dissipation occurs pre- 
dominantly in the outer 2 x 10 -5 in mass (reaching down 
to ~ 2 x 10 -4 ), or to 10% if it occurs predominantly in 
the outer 5 x 10~ 6 (reaching down to ~ 5 x 10~ 5 ). Data 
for Earth show that 1% of the absorbed solar radiation 
is converted to kinetic energy and dissipated in the at- 
mosphere, and 1% is plausible for Pegasi planets too. To 
alter the evolution, the energy need be deposited only a 
few scale heights below the altitude where it is created, 
lending plausibility to the idea. 

The presence of energy dissipation may be quantified 
in the future when several Pegasi planets have been char- 
acterized. With several ground programs (STARE, VUL- 
CAN), accepted space missions (COROT, MONS, MOST) 
and proposed ones (KEPLER, EDDINGTON) aiming at 
detecting photometric transits of Pegasi planets, there is 
indeed a good chance that enough statistical information 



on the mass radius relationship of Pegasi planets can be 
gathered. 

An unfortunate consequence of this study is that the 
possibility to determine the planets' compositions solely 
from their mass, radius and orbital characteristics seems 
to be postponed to a more distant future. On the other 
hand, we should rejoice over the perspective of better un- 
derstanding of irradiated atmospheres and tidal dissipa- 
tion. As usual, progress will mainly occur through obser- 
vations and the direct characterization of Pegasi planets. 
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